Aging Correlation Functions for Blinking Nano-Crystals, and Other On - Off 

Stochastic Processes 
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Following recent experiments on power law blinking behavior of single nano-crystals, we calculate 
two-time intensity correlation functions (I(t)I(t + 1')) for these systems. We use a simple two state 
(on and off) stochastic model to describe the dynamics. We classify possible behaviors of the 
correlation function and show that aging, e.g., dependence of the correlation function on age of 
process t, is obtained for classes of the on time and off time distributions relevant to experimental 
situation. Analytical asymptotic scaling behaviors of the intensity correlation in the double time 
t and t' domain are obtained. In the scaling limit (I(t)I(t + t')) — > h(x), where four classes of 
behaviors are found: (i) finite averaged on and off times x = t! (standard behavior) (ii) on and 
off times with identical power law behaviors x = t/t' (case relevant for capped nano-crystals). (iii) 
exponential on times and power law off times x = tt (case relevant for uncapped nano-crystals). 
(iv) For defected off time distribution we also find x — t + t' . Origin of aging behavior is explained 
based on simple diffusion model. We argue that the diffusion controlled reaction A + B f± AB, 
when followed on a single particle level exhibits aging behavior. 



I. INTRODUCTION 

The fluorescence emission of single colloidal nano- 
crystals (NC), e.g. CdSe quantum dots, exhibits interest- 
ing intermittency behavior Q. Under laser illumination, 
single NCs blink: at random times the NC will turn from 
state on in which many photons are emitted, to state off 
in which the NC is turned off. One method to character- 
ize blinking quantum dots is based on the distribution of 
on and off times. According to the theory of Efros and 
Rosen jjj, these on and off times, correspond to a neu- 
tral and ionized NC respectively. Thus statistics of on 
and off times teaches us on ionization events on the level 
of a sing le NC. Surprisingly, @,H distributions of on and 
off times exhibit power law statistics. For capped NCs 
the probability density function (PDF) of on time decays 
likeV+W ~ t-^ a +, while for off times ip-(t) ^i" 1 ""-, 
where in many cases a+ and a_ are close to 0.5 0. 

Statistical behavior of single emitting NCs, and more 
generally single molecules 0| or atoms 0,0, [Hl> ls usu_ 
all y ch aracterized based on intensity correlation functions 
[rLTTiL ItI ] . The calculation of intensity correlation func- 
tions, and the related Mandel Q parameter, for single 
molecule spectroscopy is a subject of intense theoretical 
research El El El El El El El HI El E2 (see H for 
review). Experiments on single NCs show how the cor- 
relation function method yields dynamical information 
over time scale from nano-second to tens of seconds Q. 
The correlation function of single NCs exhibits a non- 
ergodic behavior, as such these systems exhibit behavior 
very different than other single emitting objects. 

The goal of this paper is to calculate the averaged in- 
tensity correlation function for the emitting NCs. For 
this aim we use a simple two state stochastic model. The 
motivation for the calculation is twofold. First, the av- 
eraged correlation function exhibits interesting aging be- 
havior, as we will demonstrate. This aging behavior is a 
signal of non-ergodicity. Secondly, to obtain understand- 



ing of non-ergodic properties of the correlation function, 
one must first understand how the averaged correlation 
function behaves. In a future publication we will discuss 
the non-ergodic behavior of the NC correlation function, 
namely the question of the distribution of correlation 
functions obtained from single trajectory measurements. 

Aging in our context means that the (non-normalized) 
intensity correlation function 



C(t,t') = {I(t)I(t + t')), 



and the normalized correlation 



\t,t') 



C(t,t') 



(I(t)){I(t+t')) {I(t)){I(t+t')) 



(1) 



(2) 



depend on the the age of the process t even in the limit 
of long times. Here I(t) is the fluctuating stream of 
photons emitted from the NC (units counts per second). 
In the ergodic phase (i.e. when both the mean on and 
off times are finite) stationarity is reached meaning that 
C(t,f) -> C(f) when t -> oo and similarly for g^(t,t'). 
The average in (QJ is over many single NC intensity tra- 
jectories. 

Previously, Jung, Barkai and Silbey [24[ showed the re- 
lation of the problem to the Levy walk model [23|. The 
approach in j2J| is based on the calculation of MandePs Q 
parameter and does not consider the aging properties of 
the NCs. Verberk and Orrit [11] considered the problem 
of the correlation function for blinking NC, however they 
assume that the mean on and the mean off times are 
finite, while the experiments show an infinite off and on 
times (for capped NCs). To overcome this problem Ver- 
berk and Orrit introduce cutoffs on the on and off times. 
The results of Verberk and Orrit are different than ours: 
they do not exhibit aging and they are meant to describe 
the correlation function of a single trajectory (however 
the ergodic problem was not considered). Brokmann et 
al. 0| have measured aging behavior of a number of NCs. 
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They concentrate on the measurement of the persistence 
probability (see details below) while this work is devoted 
to the investigation of the intensity correlation function. 

We note that concepts of statistical aging and persis- 
tence, used in this manuscript, were introduced previ- 
ously in the context of the trap model an d g lassy dy- 
namics by Bouchaud and co-workers 

00 El & Sta " 
tistical aging is found in continuous time random walks 

pul E3|, and in deterministic dynamics of low dimen- 
sional chaotic systems j^. Aging in complex dynamical 
systems, for example super-cooled liquids or glasses is a 
topic of much research [2J|. In contrast we will later show 
that aging in NCs may be a result of very simple physical 
processes (e.g. normal diffusion). Thus we expect aging 
and non-ergodic behavior to be important in other single 
molecule systems. 

In the context of fractal renewal theory, Godreche 
and Luck [35] have considered the problem of the av- 
eraged correlation, however, in the language of single NC 
spectroscopy, they assume that statistical properties of 
the on time are identical to the statistical properties 
of off times, i.e. = ^-(t). Here we use meth- 
ods developed in [^H to the case relevant to experiments 
ip+(t) ^ ip-(t). We also obtain the aging correlation 
function in the scaling limit in the time domain. 

This paper is organized as follows. In Sec. [H]the math- 
ematical model is presented and the physical meaning of 
on and off times distributions is discussed. A brief dis- 
cussion of ensemble average and time average correlation 
function is given. In Sec. IIIII statistical properties of the 
stochastic process are considered, e.g. average number 
of jumps etc. In Sec. IIVI the distribution of the forward 
recurrence time is calculated, the latter is important for 
the calculation of the aging correlation function. In Sec. 
IVl we calculate probability of number of transitions be- 
tween t and t + t', with which the mean intensity (Sec. 
IVIJI and the aging correlation function (Sec. IVIIJl are 
obtained. Sec. I Villi is a summary. 



II. STOCHASTIC MODEL, AND DEFINITIONS 

The random process considered in this manuscript, is 
schematically depicted in Fig. The intensity I(t) 

jumps between two states I(t) — +1 and I(t) — 0. At 
start of the measurement t = the NC is in state on: 
1(0) = 1. The sojourn time 7$ is an off time if i is even, 
it is an on time if i is odd (see Fig. QJ- The times Tj 
for odd [even] i, are drawn at random from the probabil- 
ity density function (PDF) respectively. 
These sojourn times are mutually independent, identi- 
cally distributed random variables. Times t% are cumu- 
lative times from the process starting point at time zero 
till the end of the i'th transition. Time t on Fig. ^is the 
time of observation. 

We denote the Laplace transform of ip± (t) using 



V»±(s) 



1>±(t)t 



£ dt. 



(3) 




', «2 <3 



N t 'N+1 



Figure 1: Schematic temporal evolution of the dichotomous 
intensity process. 



We will classify behaviors of observables of interest using 
the small s expansion of ip± (s). We will consider: 
(i) Case 1 PDFs with finite mean on and off times, 
whose Laplace transform in the limit s — > satisfies: 



i>±{s) = l — ST± + 



(4) 



Here t + (t_) is the average on (off) time. For example 
exponentially distributed on and off times, 



4>±(s) 



I + ST± 



(5) 



belong to this class of PDFs. 

(ii) Case 2 PDFs with infinite mean on and off times, 
namely PDFs with power law behavior satisfying 



ip± oc t 



a- < a + < 1, 



(6) 



in the limit of long times. The small s behavior of these 
family of functions satisfies 



^±(s) = 1 - A±s a± 



(7) 



where A± are parameters which have units of time" . We 
will also consider cases where on times have finite mean 
(a + = 1) while the off mean time diverges (a_ < 1) 
since this situation describes behavior of uncapped QD 
0. 

(iii) Case 3 PDFs with infinite mean with a+ = cn_ = a 



l-A+s a 



(8) 



Note that Brokmann et al. jy| report that for CdSe dots, 
Q+ = 0.58 ± 0.17, and a_ = 0.48 ± 0.15, hence within 
error of measurement, a ~ 0.5. 

(iv) In Sec. IVII PI we will briefly consider the behavior 
of the correlation function for defected ip- (t) . 
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A. Physical Meaning of ip±(t) 

As mentioned in Introduction, and following Ref. Q, 
we assume that a charged (neutral) uncapped NC is 
in state off (on), respectively. Physically, for charged 
NC, Auger non-radiative decay time of a laser excited 
electron-hole pair, is much faster than the radiative time 
of the electron-hole pair [2fJ. Hence a charged NC is in 
off state. The physical mechanism responsible for the 
power law blinking (i.e. charging) behavior of NCs is 
still unclear. Models based on trapping of charge car- 
riers in the vicinity of the NC, and fluctuating barrier 
concepts were suggested in 0, 0, E|- Here we will em- 
phasize an alternative simple picture based on diffusion 
concepts. Before further experiments are performed, it is 
impossible to say if the simple picture we consider here 
works better or worse than other approaches. 

We note that the simplest diffusion controlled chemi- 
cal reaction A + B f± AB, where A is fixed in space, can 
be used to explain some of the observed behavior on the 
uncapped NCs. As mentioned the latter exhibit exponen- 
tial distribution of on times and power law distribution 
of off times. The on times follow standard exponential 
kinetics corresponding to an ionization of a neutral NC 
(denoted as AB). A model for this exponential behavior 
was given already in [2j . Once the NC is ionized (A + B 
state) we assume the ejected charge carrier exhibits a 
random walk on the surface of the NC or in the bulk. 
This part of the problem is similar to Onsager's classical 
problem of an ion pair escaping neutralization (see e.g., 
|3flL ErTp. The survival probability in the off state for 
time t, S-(t) is related to the off time distribution via 
S-(t) = l-J^.(r)dr, or 



M*) = - 



dS-(t) 
alt 



(9) 



It is well known that in three dimensions survival prob- 
ability decays like t -1 / 2 , the exponent 1/2 is close to the 
exponent measured in the experiments. In infinite do- 
main the decay is not to zero, but the 1/2 appears in 
many situations, for finite and infinite systems, in com- 
pletely and partially diffusion controlled recombination, 
in different dimensions, and can govern the leading be- 
havior of the survival probability for orders of magnitude 
in time ^ilElG3]- I" this picture the exponent 1/2 does 
not depend on temperature, similar to what is observed 
in experiment. We note that it is possible that instead of 
the charge carrier executing the random walk, diffusing 
lattice defects which serve as a trap for charge carrier are 
responsible for the blinking behavior of the NCs. 

One of the possible physical pictures explaining blink- 
ing of capped NCs can be based on diffusion process, 
using a variation of a three state model of 0] . As men- 
tioned in Introduction, for this case power law distri- 
bution of on and off times are observed. In particular, 
neutral capped NC will correspond to state on (as for un- 
capped NCs). However, capped NC can remain on even 
in the ionized state. We assume that the ionized capped 



NC can be found in two states: (i) the charge remain- 
ing in the NC can be found in center of NC (possibly a 
de-localized state), (ii) charge remaining in the NC can 
be trapped in vicinity of capping. For case (i) the NC 
will be in state off, for case (ii) the NC will be in state 
on. The main idea is that the rate of Auger nonradiative 
recombination Q of consecutively formed electron-hole 
pairs will drop for case (ii) but not for case (i). We note 
that capping may increase effective radius of the NC, or 
provide trapping sites for the hole (e.g., recent studies by 
Lifshitz et al. [43| demonstrate that coating of NCs cre- 
ates trapping sites in the interface). Thus the off times 
occur when the NC is ionized and the hole is close to the 
center, these off times are slaved to the diffusion of the 
electron. While on times occur for both a neutral NC 
and for charged NC with the charge in vicinity of cap- 
ping, the latter on times are slaved to the diffusion of 
the electron. In the case of power law off time statistics 
this model predicts same power law exponent for the on 
times, because both of them are governed by the return 
time of the ejected electron. 

The main point we would like to emphasize is that 
several simple mechanisms might be responsible for the 
power law statistics, and hence aging correlation func- 
tions in single molecule experiments may turn out to be 
wide spread. Beyond single molecul e sp ectroscopy we 
note that certain single ion channels |4lL E^L , deter- 
ministic diffusion in chaotic systems j4fl| , the sign of mag- 
netization of spin systems at criticality j^|, all exhibit 
intermittency behavior, and the correlation function we 
obtain here might be useful also in other fields. Hence 
we don't restrict our attention to the exponent 1/2, as 
there are indications for other values of a between and 
1, and the analysis hardly changes. 



B. Definition of Correlation Functions 

Since the process under investigation is non-ergodic, 
and since measurements are made on a single molecule 
level, care must be taken in the definition of averages. 
From a single trajectory (ST) of I(t), recorded in a time 
interval (0, T), we may construct the time average corre- 
lation function 



C-'' i(t + t')i(t)dt 



C ST (T,t') = 



T-t' 



(10) 



On the other hand we may generate many intensity tra- 
jectories one at a time, to obtain C{t,t') and g( 2 >(t,t'). 
We call C(t,t') single particle averaged correlation func- 
tion. For non-ergodic processes Cst{T, t') ^ C(t, t') even 
in the limit of large t and T. Moreover for non-ergodic 
processes, even in the limit of T — > 00, CsT{T,t') is a 
random function which varies from one sample of I(t) to 
another. 

We stress that the single particle averaged correlation 
function is not the correlation function obtained from 
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measurement of an ensemble of particles. To see this 
consider the intensity of N blinking NCs, 



'A' 



N 



(ii) 



and j is an index of the particle number. The corre- 
sponding normalized correlation function is 



G (2) (t f/) = (iN(t)I N (t + t')) 

(i N (t))(i N (t+t')y 



(12) 



If the blinking behavior of individual NCs is independent, 
but statistically identical, 

(2) , _ N(I(t)I(t + 0) + N(N - l)(I(t))(I(t + tj) 
[ ' } N*(I(t))(I{t + t>)) 

(13) 

Clearly for N = 1 we obtain the correlation function 
C(t,t') in Eq. Q, while for N — > oo the intensity does 
not fluctuate at all (as well known) . Hence while we con- 
sider here average over an ensemble of trajectories, we 
are not reconstructing the correlation function obtained 
from a measurement of a large number of NCs. Our the- 
ory is valid only for single molecule measurements, and 
the aging behavior of the normalized correlation function 
cannot be obtained from macroscopic measurement. 



III. NUMBER OF JUMP EVENTS BETWEEN 
AND t 

In this Section we investigate basic statistical proper- 
ties of the on-off process. 

The probability of n transitions (either off — ► on or 
on — > off ) between times and t is 

P t (n) = P(0, t, n) = (9 (t n < t < t n+1 )) (14) 

where 9(t n < t < i„+i) is 1 if the event in the parenthesis 
occurs; otherwise it is zero. Laplace transforming Eq. 
Ill4t with respect to t yields 



t„+i 



die 



st\ _ / a -st n 



') = (e 



-ST„ + i 



)• (15) 



A simple calculation using t n = J2i=i T * yields 



Ps (n) = 



t/j + (s) (s) 



i/2 



t[} + (s) ' ip- (s) 



2 1— 



n even 



n odd 



(16) 

To derive Eq. itTfif we used the initial condition that the 
state of the process at t — is +. Eq. 1)1611 satisfies the 
normalization condition Y^=o Ps{ n ) = V s - 



A. Mean number of renewals (n) 

Using Eq. ill fit the mean number of transitions is in 
Laplace t — > s space 



<Ms) i + ^-00 



l-^_(s)V+(s) 



(17) 



Using Eqs. QCHHJ, the small s expansion of Eq. ifTTIl . 

and then inverting to time domain we get the long time 
behavior 



(n(t)> 



21 



T + +T- 

2t a - 

A_r(i+ Q _) 

2t°_ 

I (A + +A_)r(l+a) 



case 1 

case 2 (18) 
case 3. 



B. Asymptotes of Pt (n) 

For narrow PDFs, i.e. case 1, and for long times we 
obtain from Eqs. 10] 11 lit 



Pt{n) 



{ 



n even 



< 



r +) X ( n t 
TFZJ d {2 - <r + ) + (r_) 

^(^-T^fe) n odd 



(19) 

To obtain this result we used the small s expansion of Eq. 
(tl6l and then a simple Laplace inversion. We neglected 
the fluctuations in this treatment, the latter are expected 
to be Gaussian in the long time limit. 

For broad PDFs satisfying a + = a- = a, case 3, we 
find 



P s (n) ~ A±s a - 1 e- k{A++A - )sa 



(20) 



where k — n/2 and A± = A + for n even, while k = 
(n + l)/2 and A± = A- for n odd. Inverting to the time 
domain we find 



Pt(n) 



A, 



t 



(A+k + A-k) 1 /^ 1 a 



t 



(A+k + A-k) 



l/a 



(21) 

where l a (t) is the one sided Levy stable PDF whose 
Laplace pair is exp(— s Q ). 

For case 2, with a_ < aq_, we get 



P s (n) exp(-A_/fcs Q -) 



(22) 



for k = n/2 and n even. The probability of finding an 
odd n in this limit is zero. This is expected since the off 
times are much longer than the on times, in statistical 
sense. Thus for long times we have 



Pt(n) 



A. 



t 



a- {A-k) 1 / a -+ 1 
and n is even. 



t 



(A_k) 



(23) 
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IV. FORWARD RECURRENCE TIME 

The time E — tN+i—t is called the forward recurrence 
time. The times (see Fig. t^+i and t/v are defined in 
such a way that tjsr <t < tjv+i, hence also N is a random 
variable. Let ft(E) be the probability density function 
of the random variable E. The subscript t in ft{E) indi- 
cates that t is a parameter, while E is a random variable. 
Generally the PDF of E depends on how old the process 
is, namely on t. A process is said to exhibit statistical 
aging if even in the limit of t — > oo, ft(E) depends on 
t. The PDF ft(E) is important for the calculation of the 
aging correlation function. 

We consider the joint PDF 

i t (E, N) = {5{E- t N+1 +t)6{t N <t< tN+i)). (24) 

Later we will sum over N to obtain ft(E). We consider 
the double Laplace transform C of Eq. ll24T) with t —> s 
and E — > u 

C t . E {it{E,N)} = 



At 



dEe 



-st-uE 



6(E-t N+1 +t)) 



— (e _u * JV + 1 / e ~( s - u )t 



di) 



— (s— u)tisr+i _ p —(s—u)t[ 

— ^"*JV+l. 

u — s 

For even N we use the averages 



(25) 



and find 



f,(u,JV) = 



(e- u[tN + 1 ' tN) ) =i)+ («). 



JV/2 

"0+ (s) "0- («) ^+ (s) - V>+ («) 



w — s 

In similar way we obtain for N odd 



(26) 



U{u,N) 



0+ ( s ) 0- ( s ) i>-(s) -i>-(u) 4>+(s) 



(27) 



Note that { s (u,N) is the double Laplace transform of 
^(-E 1 , AT), while t/>±(s) and 4>±{u) are single Laplace trans- 
forms. Summing f s (u,N) over TV we obtain the double 
Laplace transform of /( (I?) 



/,(«) = 







-0_(s) - ^_ (ti) 


(u - s) 


1 





(28) 



For ^+(s) = 0_(s) Eq. reduces to Eq. 6.2 in [35| . 



A. Limiting cases for ft(E) 

We now analyze the long time t — > oo behavior of 
ft(E). In this case we expect that an equilibrium PDF 
for ft{E) will emerge. This equilibrium is related to sta- 
tionarity, ergodicity, and aging as we will show. 

We consider narrow distributions, i.e. case 1 first. Tak- 
ing the limit s — > of Eq. H28JI . corresponding to t — > oo 
and find 



fs{u) 



1 — ip + (u) 1 — ?/>_ (u) 



SU (t+ + T_ ) ,SU (t + + T_ ) 



(29) 



The Laplace s — y t and u — > 15 inversion of this equation 
is immediate 



(30) 



a behavior which is valid in the limit of long time t (and 
independent of it). Note that the first (second) term 
on the right hand side of the equation, corresponds to 
trajectories with even (odd) number of steps. One can 
show that in the limit of long times probability of finding 
the process in state ± is 



lim P±(t) = 

f— >O0 T+ + T_ 



(31) 



as might be expected. In the special case of tp- = ip+ 
we obtain a well known equation which has several 
applications in theory of random walks e.g. [131 ■ The 
important point to notice is that in the limit of large time 
i, and when average times t± are finite, an equilibrium 
is obtained which does not depend on t. 

We now consider broad distributions, with diverging 
averaged on and off times, case 2. In the limit of small 
s and small u. with their ratio finite 



/.(«) 



(32) 



The investigation of this equation yields the long time t 
behavior 



f t (E) ~ d t (E), 



(33) 
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where dt(E) is Dynkin's function 



d t (E) 



sin(7ra_) 



E a -(t + E)' 



(34) 



From Eqs. 1321331 we learn that unlike case 1, the PDF of 
E depends on time t even in the long time limit. Eq. 13411 
was obtained by Dynkin pV?! l3^ | as a limit theorem for 
renewal processes with a single waiting time PDF. Here 
we showed that for a two state process the details on a + 
and A + are not important in the long time limit. This 
is expected, the off (i.e. minus) times are much longer 
than the on (i.e. plus) times in statistical sense, and 
hence our results in the long time limit are not sensitive 
to the details of i>+(t). In the same spirit it can be shown 
that in the limit of long time t, and with probability one, 
the process is found in state minus. 

Finally for case 3 where a+ = a_ < 1, we find that 
Eqs. 1321331 are still valid. However now probability of 
finding the process in state ± is given by 



P± = lim P±(t) 



A, 



A+ + A- 



(35) 



B. Joint PDFs for Forward Recurrence time 

It will turn out important to define the joint PDFs of 
time E provided that process is in state plus or state mi- 
nus at time t. We denote these PDFs with ft(E, ±) and 
the corresponding double Laplace transform f s (u,±). 
Since the start of process is state + at time t = 0, we 
get using Eq. ll2"6T) 



fs(u,+) 



n=o,n even 



f s (u,N) 



(u-s) 1 - ip + (s)ip-(s) 
and using Eq. ifSTjl 



f s (u,-)= Uu>n) 

n=o,n odd 



(36) 



*!>+(*) 



(u - s) 



l-V+(i#-(s) 



Note that 



f s (u) = f s (u,-) + f s (u,+). 



(37) 



(38) 



The probability of finding the particle in state + when 
t — > oo is 



/-■oc 

lim P+ = lim / f t (E,+)dE, 

t^OO t — J-OO Jq 



(39) 



provided that the limit exists. For example for case 1 it 
is easy obtain from Eq. l(39jl the result in Eq. (plj. 

The limiting PDFs f t {E,±) are obtained in double 
Laplace space by considering the small s (and small u 
for cases 2 and 3) limit. They are 



fs(u,±) 



ms(t + +t_) 

A± u a —s a 
A + +A_ (u — s)s° 



case 1 



case 3. 



(40) 



For case 2 we find in this limit f s (u, +) = 0, i.e. proba- 
bility of finding the particle in state + is zero, and 



fs(u, -) 



(u — s)s a 



(41) 



The double inverse Laplace transform of this equation is 
given in Eq. (l3"4l) . 



V. NUMBER OF RENEWALS BETWEEN TWO 
TIMES 

We now calculate P(t,t + t' , n) the probability of num- 
ber of renewals n between time t and time t + t'. Ob- 
viously the process is generally not stationary and the 
information on P(0,t',n), obtained in Sec. IIII1 is not 
sufficient for the determination of P(t, t + t',n). We now 
classify the trajectories according to the state of the pro- 
cess (i.e., + or — ) at times t and t + t' . It will turn 
out that the intensity trajectories, when the process is in 
state + at time t and state + at time t + i', are those 
which are important for the calculation of the correlation 
function. 

The probability of not making a jump in time interval 
t,t + t', when the process is in state + at time t and state 
+ at time t + t' is 



P++(t,t 



t',0) = J f t (E,+)dE. 



(42) 



The probability of finding n > transition events in time 
interval t,t + t' ', when state of process at time t is + and 
state of process is t + t' is also + 



C t 'P ++ (t,t + t',n) 



ft(u,+)i>-(u) 4>+(u)->Jj-(u) 



l-^+(«) 



(43) 



where n > is even, and u is the Laplace conjugate of 
t' . Note that Eq. jHH depends on t through f t (u,+). 

The other combinations, e.g. P-\ (t,t + t' ,n), are given 

in AppendixEI as we U as P(t> t + t', n). 



VI. MEAN INTENSITY OF ON-OFF PROCESS 

The averaged intensity (I(t)) for the process switching 
between 1 and and starting at 1 is now considered. In 
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Laplace t — > s space it is easy to show that 

i-iM«) i 



l-V+(s)V-(s) 



(44) 



One method to obtain this equation is to note that 
(J(t)) = Prob[J(t) = 1], hence 



I(s) 



f a (E,+)dE, 



(45) 



and therefore the u — * limit of Eq. ERfo yields Eq. (14*41 . 

The Laplace s — »■ t inversion of Eq. (33} yields the 
mean intensity (-?(£)). Using small s expansions of Eq. 
144f , we find in the limit of long times 



(A*)> 



T++T- 

A + t"-~' 



A_r(l+a_-a + ) 



case 1 
case 2 
case 3. 



(46) 



If the on times are exponential, as in Eq. (J"} then 



7( S ) 



1 + ST + - 1p_ (s) 



(47) 



This case corresponds to the behavior of the uncapped 
NCs. The expression in Eq. (3ZJj and more generally, 
the case a_ < a + = 1 leads for long time t to 



(Ht)} 



T+t c 



^_r(a_) 



(48) 



For exponential on and off time distributions Eq. (§}, 
we obtain the exact solution 



r_ exp 



l 



(49) 



Remark For the case a+ < a_ < 1, corresponding to 
a situation where on times are in statistical sense much 
longer then off times, (I(t)) ~ 1, 



correlation function. Summing Eq. 14-tll over even n > 0, 
and using Eq. (4*2)) for n = 0, we find 



C(t,u) = 



A(u = 0,+)-/ t (u,+) 



+ /t («,+)- 






-(«) 


1-^4 


(«)' 




1 — lj}_{u)lj}. 





(50) 



where u is the Laplace conjugate of t'. We see that the 
correlation function generally depends on time t. 



A. Case 1 

For case 1 with finite r + and r_, and in the limit of 
long times t, we find 

lim C(t, u) = 

t — >oo 



1 - ± , 5- (51) 



T + 1i 



This result was obtained by Verberk and Orrit 0] and 
it is seen that the correlation function depends asymp- 
totically only on t' (since u is Laplace pair of t'). Namely, 
when average on and off times are finite the system does 
not exhibit aging. If both tjj + (t) and ip- (t) are exponen- 
tial then the exact result is 



C(t,f) = 



r_ exp —tl 



+ T-, 



T_ exp[-t'(^ + ^jj+r + 



and C(t,t') becomes independent of t exponentially fast 
as t grows. 



B. Case 2 



VII. AGING CORRELATION FUNCTION OF 
ON- OFF PROCESS 

We are now able to calculate the correlation function 
C(t,t') = (I(t)I(t + t')). We consider the process I(t) 
as jumping between state on with I{t) = 4-1 and state 
°ff I{t) = 0- The symmetric case where I(t) jumps be- 
tween the states I(t) = —1 or I(t) = —1, is discussed in 
Appendix [B| We assume that sojourn times in state on 
(off) are described by ip+(t) (ip-(t)), respectively. Con- 
tributions to the correlation function arise only from tra- 
jectories with I(t) = 1 and I(t + t') = 1, meaning that 
only the 4-+ trajectories, in Eq. 14-'tt contribute to the 



We consider case 2, however limit our discussion to 
the case a + = 1 and a_ < 1. As mentioned this case 
corresponds to uncapped NCs where on times are expo- 
nentially distributed, while off times are described by 
power law statistics. Using the exact solution Eq. (50)1 
we find asymptotically, when both both t and t' are large: 

Unlike case 1 the correlation function approaches zero 
when t — > oo, since when t is large we expect to find 
the process in state off. Using Eq. (3BI> the asymptotic 
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behavior of the normalized correlation function Eq. <|2J 
is 

5 (2) (i,t')-(l + ^) 1_a ". (53) 

We see that the correlation functions Eqs. l)52l l53l) ex- 
hibit aging, since they depend on the age of the process 
t. 

Considering the asymptotic behavior of C(t,t') for 
large t, but small £', yields in the limit of s <C u, s <C 

C(t,u) pa 



o' 




1 



l — lj) + {u) 1 — i}-{u) 



(54) 

This equation is similar to Eq. 1)5111 . especially if we 
notice that the "effective mean" time of state off until 
total time t scales as A-t 1- "- . Despite assumptions of 
s <C u, s <C uip(u) in the derivation of Eq. 115411 . it also 
reproduces the result Eq. H52I1 and hence is applicable 
for any u (and thus t 1 ) as long as t is large enough. 

For the special case, where on times are exponentially 
distributed, the correlation function C is a product of 
two identical expressions for all t and t 1 : 



C{s,u) 



1 + ST + — (s) 1 + UT + — 1p_ (u) ' 



(55) 



where s (u) is the Laplace conjugate of t (f) respectively. 
Comparing to Eq. (T7)l we obtain 



C(t,t') = {I(t)){I(t% 



and for the normalized correlation function 



(Ht')) 

(i(t + t')y 



(56) 



(57) 



Eqs. 1(571 l5fijl are important since they show that mea- 
surement of mean intensity (I(t)) yields the correlation 
functions, for this case. While our derivation of Eqs. 1)571 
l56| is based on the assumption of exponential on times, 
it is valid more generally for any tp+(t) with finite mo- 
ments, in the asymptotic limit of large t and t' . To see 
this note that Eqs. ^EH) yield C(t,t') - (I{t))(I(t'). 

In Fig. [21 we compare the asymptotic result ll52|l with 
exact numerical double Laplace inversion of the corre- 
lation function. We use exponential PDF of on times 
ip+{s) = 1/(1 + s), and power law distributed off times: 
i>-(s) = ip-(s) = 1/(1 + s 0A ) corresponding to a_ = 0.4. 
Convergence to asymptotic behavior is observed. 

Remark For fixed t the correlation function in Eq. 
ll52ll exhibits a (t') 01 -^ 1 decay. A (t') 01 -^ 1 decay of an in- 
tensity correlation function was reported in experiments 



Figure 2: Exact C(t,t') for Case 2: exponential on times 
and power law off times with a_ = 0.4. We use ^+(s) = 
1/(1 + s) and ip-(s) = 1/(1 + s 0,4 ) and numerically obtain 
the correlation function. For each curve in the figure we fix 
the time t. The process starts in the state on. Thick dashed 
straight line shows the asymptotic behavior Eq. 15211 . For 
short times (t' < 1 for our example) we observe the behavior 
C(t,t') ~ C(t, 0) = (I(t)), the correlation function is flat. 



of Orrit's group 0] for uncapped NCs (for that case 
a_ = 0.65 ± 0.2). However, the measured correlation 
function is a time averaged correlation function Eq. 11 Oil 
obtained from a single trajectory. In that case the cor- 
relation function is independent of i, and hence no com- 
parison between theory and experiment can be made yet. 



Case 3 



We now consider case 3, and for long t and t' we find 



„, ^ sm7ra 
C(t,t') = P+- P+P- B 



7T \l+t/t' 
with P± given by Eq. II85H and where 



1 — a, a 



(58) 



B(z; a, b) 



is the incomplete beta function. The behavior in this 
limit does not depend on the detailed shape of the PDFs 
of the on and off times, besides the parameters A + /A- 
and a (see also Appendix Q- We note that both terms 
of Eq. il5fill contribute to Eq. 1(58)1 . The appearance 
of the incomplete beta function in Eq. 1158)1 is related 
to the concept of persistence. The probability of not 
switching from state on to state off in a time interval 
(t, t + t'), assuming the process is in state on at time t, 
is called the persistence probability. In the scaling limit 
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Figure 3: Exact C(t,t') for case 3, when both on and off times 
are power law distributed with a — 0.4. We use ip±{s) = 
1/(1 + s 0A ) for different times t increasing from the topmost 
to the lowermost curves. The dots on the left and on the right 
show C(t,0) = (/(*)) and C(t,ao) = (/(*)) /2 respectively. 
The process starts in the state on. 



this probability is found using Eqs. 1-141 14()|) : 

sin(Tra) f°° t a 



E a (t + E) 



dE 



1 



sm 7ra 



-D 



1 + t/t' ' 



1 



(59) 



The persistence implies that long time intervals in which 
the process does not jump between states on and off, 
control the asymptotic behavior of the correlation func- 
tion. The factor P + , which is controlled by the ampli- 
tude ratio A+/A-, determines the expected short and 
long time t' behaviors of the correlation function, namely 
C(oo,0) = lim t ^ oo (7(t)J(t + 0)) = P+ and C(oo,oo) = 
lim t ^ 00 (7(t)J(t + oo)) = {P+f. With slightly more de- 
tails the two limiting behaviors are: 



P + 



7«! 



C{t,t>) 



sin(7ra) 

TTOi 



(60) 



*e > I- 



Using Eq. 14fill the normalized intensity correlation func- 
tion is <? (2) (M') ~C(t,t')/(P+) 2 . 

In Fig. 03 we compare the asymptotic result II58J1 with 
exact numerical double Laplace inversion of the correla- 
tion function for PDFs $+(s) = ^_(s) = 1/(1 + s 0A ). 
Convergence to Eq. ll58l) is seen. 

Remark For small t'/t we get flat correlation func- 
tions. Flat correlation functions were observed by Da- 
lian's group Q for capped NCs. However, the mea- 
sured correlation function is a single trajectory correla- 
tion function Eq. iflOl) , and hence no comparison between 
theory and experiment can be made yet. 



D. Defected off time distribution 

As mentioned in the Introduction, an off state of un- 
capped NC corresponds to an ionized NC. Assume that 
the transition from state on to state off occurs when a 
charge carrier is ejected into the vicinity of the NC, and 
then starts to move diffusively in the bulk. If the diffusion 
process takes place in three dimensions, there is a finite 
probability that the charge carrier will not return to the 
NC. In that case the NC remains in state off forever. 

Such a situation can be modeled based on defected 
distribution of off times. In this case we have a non- 
normalized PDF of off times 



ip^(t)dt = Z < 1, 



(61) 



the small s expansion of the Laplace transform of i/j- (t) is 
^>(s) ~ Z — A-s a - . Z is the probability of charge carrier 
to return; this probability was the subject of extensive 
investigation in the context of first passage time problems 

na 

For large t the mean intensity is 



(/(*)> 



A + t~° + 

(i-z)r(i-« + ) ' 

A + A_t-( 1+ °-> 
(l-Z)2|r(-a_)|> 



a± < 1 



ev_ < a i = 1. 



(62) 



Note that here a_ can be smaller, larger or equal to a+ 
when a+ < 1. 

Using Eq. (IBTH) we obtain asymptotically, for both t 
and t' large, 



C{t,t') 



(i-z)r(i- Q+ ) L ) ' 



a± < 1 



(i-zm-a-) ; (tt')- {1+a -\ a_<a + = l. 

(63) 

Using Eq. Ki2ll we can relate the intensity correlation 
function with the mean intensity, 



C(M') 



(I(t + f)), a±<l 
(I(t))(I(t')), a_<a+ = l. 



(64) 



This result shows that for a + < 1 (I(t)I(t + 1')) ~ 
(I(t + t')), independent of the value of I(t). Asymptotic 
validity of this relation can be explained by noticing that 
the non-zero contributions to (I(t)I(t + 1')) come only 
when both I(t) and I(t + t') are equal to 1. However, at 
long time t, after jumping off there is a negligible prob- 
ability of being on again for a cumulative time duration 
comparable with the total time (note that (I(t)) scales as 
the probability of making no transition off, i.e., of persis- 
tence). Hence, the nonzero contributions to (I(t + t')) 
are mainly those staying on from time t, so that al- 
most certainly, if I(t + 1') — 1 then also I(t) = 1 and 
(I(t)I(t + t>))~ (I(t + t>)). 
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The above argument does not hold in the case of a+ — 
1, because now the on times are very short and there 
is no possibility of staying on persistently for long times. 
Accordingly, the decay of (I(t)} and C(t, t') is much faster 
here. The leading contributions to (I(t)) and C(t,t') in 
the case of a + < 1 disappear as a+ f 1, due to the pole 
of the T-function. 



VIII. SUMMARY 

We demonstrated the dependence of the two-time cor- 
relation function C(t,t') = (I(t)I(t + t')) on the times 
t and t'. This is in a full contrast to the well-known 
convergence of the correlation function to the station- 
ary limit which is independent of t. Such a convergence 
is found when the average on and off times are finite 
(as shown above for exponential on and off distribu- 
tions). When these times diverge non-stationary behav- 
ior is found. The non-vanishing i-dependence of the cor- 
relation function C(t,t') is known as aging. 

We obtain different modes of aging, yielding depen- 
dence of C(t,t') on the ratio t/t', product tt' and the 
sum t + 1'. 

(i) For PDF of on times having finite mean and power 
law distributed off times with infinite mean, the corre- 
lation function asymptotically splits into a product of 
two identical functions, one of t and the other of t' (see 
Eq. jlH ), leading to tt' dependence Eq. 15211 . This case 
corresponds to the behavior of the uncapped NCs. 

(ii) When both on and off times are described by broad 
distributions, with identical exponents a + = a_, the 
correlation function depends on the ratio t/t', Eq. 115811 . 
This case corresponds to the capped NCs (within the 
error of measurement). 

(iii) For defected off times and a+ < 1, we find that the 
correlation function depends on t + t' , Eq. . 

(iv) Finally, for stochastic processes with finite on and off 
times, we recover known behavior, where the correlation 
function in the scaling limit depends only on t', Eq. llolt . 

In different regimes, the correlation function exhibits 
either a strong sensitivity on the details of the stochas- 
tic process (i.e., on ip±(t)), or certain universal features 
which are now discussed. We also found relations be- 
tween the correlation function and mean intensity, for 
several cases. 

(i) For PDF of on times having finite mean and power 
law distributed off times with infinite mean, the correla- 
tion function is related to the mean intensity according 
to C(t,t') ~ (I(t))(I(t')) Eq. ldJ. For short times t' the 
correlation function depends on the details of tp± (t) , Eq. 

(S3J. 

(ii) When both on and off times are described by broad 
distributions, with identical exponents a + = a_, the per- 
sistence probability governs the aging correlation func- 
tion Eq. (EBJ. This is a universal behavior in the sense 
that all ip± (t) belonging to this family, yield identical be- 



havior for the correlation function C(t,t') in the limit of 
t — > oo. 

(iii) For defected off times and a + < 1, we find that the 
correlation function C(t, t') = (I(t + t')). 

(iv) For the standard case, where both the mean on and 
mean off times are finite, the correlation function de- 
pends on the details of tp±(t). 

Simple physical explanations for the aging behavior 
were briefly discussed. Models, discussed previously in 
the literature, based on diffusion processes, or fluctuat- 
ing barrier models, or trap models, may all lead to aging 
behavior of the correlation function. Thus aging behav- 
iors in single molecule spectroscopy may have other ap- 
plications besides single nano-crystal spectroscopy. The 
dependence of the correlation function on control param- 
eters like temperature and laser intensity can be used to 
distinguish between the microscopic scenarios proposed 
here and in the literature. 
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Appendix A: NUMBER OF TRANSITIONS 
BETWEEN t AND t + t' 



The probability of finding n > transition events in 
time interval t,t + t', when the process is in state + at 
time t, and state — at time t + 1' 



ft(u,+) ip + (u)ip-(u) 



^ 1 -$_(«) 



(Al) 



n is odd. The probability of finding n = transition 
events in time interval t,t + t', when the process is in 
state — at time t, and state — at time t + t' 



P—(t,t + t',0) = j i ft(E,-)dE, 



(A2) 



with obvious notations, and for even n 
C t 'P—(t,t + t',n) = 



f t (u, -)^+{u) il) + {u)ip^{u) 
Finally for odd n 

C t >P-+(t,t + t',n) 



1 - Tp-(U) 



(A3) 
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ft(u,-) ip + (u)ip-(u) 



r- 1 - $+(u) 



(A4) 



The indexes ij = ++, ij = , ij = — (-, and ij = + — 

correspond to trajectories which start in state i at time 
t and end with state j at time t + t' . Obviously for 
P(t, t + t' , n) we have for even n > 

£ t /P(i, t+t', n) = C v P++{t, t+t', n)+C t 'P—{t, t+t', n), 

(A5) 

while for odd n 

Lt>P{t,t+t',n) = C t 'P+-(t,t+t' ,n)+£ t >P-+(t,t+t' ,n). 

(A6) 



where i* = + if i = — and i* = — if % = +. If on and 
o/f times are identically distributed, ip+(t) = we 
obtain the result given by Godreche and Luck [35] . 



Appendix C: ASYMPTOTICS OF C(t,t') FOR 
CASE 3 (a+ = a- = a < 1) 

If, in analogy to the cases 1 and 2, we wish to explore 
the behavior of the correlation function in the limit of 
large t, but for any time difference t', it is easy to obtain 
the following asymptotic result for small s from Eq. (15011 : 



Appendix B: JUMPS BETWEEN 1 AND -1 

We now consider a correlation function which slightly 
differs than the one considered in the main text. We 
assume that I(t) = +1 or I(t) = —1, the + [— ] times are 
described by ip+ (t) [ip- (t)] respectively. For this case the 
correlation function is related to P(t, t + t', n) according 
to 



C(M') = J2{-l) n P(t,t + t',n). 



(Bl) 



In Laplace t' — *■ u space, and using convolution theorem 
of Laplace transform, 

CfC{t,t') = C(t,u) - 



Mu = 0,+)-ft(u,+) | ft(u = 0,-)-f t (u,-) 



£t>P + +(t,t+t',n)+ UP—(t,t+t',n) 



n>o, even 



n>o, even 



E 



n>o, odd 



£ t/ P-+(t,t + t',n)- 2^ Ct'P + -(t,t + t',n). 

n odd 

(B2) 

The first two terms on the right hand side of equation 
(IB2I correspond to trajectories with no transitions in 
time interval t, t + t' . The ft(u = 0, ±) are the prob- 
abilities P± of finding the process in state ± at time t. 
Using Eqs. EHEIIEIEI we get 



C{t,u) 



1 -Y,i=±ft(u,i) 



Y^ftM [&.(u)(2-&(u))-l 



»=± 



(B3) 



C(s,u) 



1 



1 



i - i>-(u) 



1 - ij)_{u)ij) + {u) 



(s - u) 

J (CI) 

However, consistent with the demand t' -C t we have 
u S> s and so have to remove s in the second term in 
the brackets of Eq. IIClJl . It would be wrong to try 
and attempt inverse Laplace inversion of this term with 
respect to s. Thus, we obtain 

C{t,t' <i) wP + , 

in agreement with Eq. llfif))) . We see that the correlation 
function is virtually constant for any t' (even small) and 
for any t large enough, as long as t' -C t. This, of course, 
could be expected based on the fact that asymptotic ex- 
pression Eq. l|58j) gives the exact value of P+ for t' = 
(see also Fig. |2J. 

Eq. ilClJl can also be used to check the asymptotics 
when t' also becomes large. Using small u expansions for 
4>±(u) yields 



C(s,u) 



P+i — 



su u(s — u) 



The double inverse Laplace transform of l/u(s — u) is 
either or 1 (for positive t,t'), depending on whether 
we assume u > s or u < s (i.e., do we first perform 
inversion with respect to u or to s). The choice u > s is 
appropriate when t' <C t and vice versa, hence we recover 
asymptotic limits shown in Eq. ijfifljl . up to the leading 
order. 

To conclude, we have demonstrated that in case 3, in 
the long t limit the correlation function does not depend 
on the particular form of tp± but only on their asymp- 
totics, for any t', as given by Eq. l|58Jl . This is in contrast 
to the long t limiting behavior of C(t,t') in cases 1 and 
2, where C(t,t') does depend on the particular form of 
ij)± for short t'. 
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